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Введение. Статья посвящена аналитическим исследованиям динамической системы «ротор — опора с зазором», 
находящейся под технологической нагрузкой. Цель исследования — получить выражения для определения эк- 
вивалентных жесткостных характеристик системы. 

Материалы и методы. Рассмотрен ротор, вращающийся в упругих опорах с зазорами. Предложена динамиче- 
ская модель, позволяющая рассматривать задачу определения линейных эквивалентных жесткостных характе- 
ристик опорных узлов. Для решения задачи составлена система дифференциальных уравнений и выполнен их 
детальный анализ. 

Результаты исследования. Из полученных уравнений динамики рассматриваемой системы можно рассчитать 
статический угол отклонения цапф ротора, обусловленный действием нагрузки. Предложенные выражения для 
определения эквивалентных жесткостных характеристик свидетельствуют о том, что можно исследовать дина- 
мику ротора как на линейных упругих опорах с указанными выше параметрами. 

Проанализирована полученная система уравнений и перечислены все частные случаи применения формул пер- 
вого приближения для эквивалентных жесткостей опорных узлов ротора. 

Обсуждение и заключения. Полученные результаты позволяют исследовать многие динамические процессы на 
основе линейных дифференциальных уравнений, учитывая нелинейные свойства системы. Для тех строгальных 
машин, которые используются при работе с кожаными материалами, определение колебаний ротора в горизон- 
тальной плоскости обеспечивает качество и точность операций. 
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Введение. Технический прогресс предъявляет новые, более высокие требования к качеству продукции 
машиностроения. Главным образом речь идет о высокой надежности, долговечности машин, их продуктивно- 
сти и безопасности. Все эти параметры должны быть учтены и рассчитаны на этапе разработки. 

С точки зрения производственных процессов особое значение имеет динамическое роторное оборудо- 
вание, обеспечивающее непрерывность технологического процесса [1-3]. В [4] представлены схемы конструк- 
ций и принцип работы различных роторных машин, выполнен патентный поиск по каждому их виду. В [5] рас- 
смотрен характер зависимости эквивалентной жесткости от частоты колебаний для определенной модели, при- 
ведена зависимость критических частот от эквивалентной жесткости опор. В [6] авторы подробно рассматри- 
вают вопросы контроля вибрации, виброналадочных работ и предотвращения увеличения вибрации роторных 
машин, равновесия роторов. Кроме того, здесь перечислены источники вибрации и приводятся базовые сведе- 
ния из теории колебаний. В работах [7-9] обсуждаются проблемы динамики жесткого неуравновешенного ро- 
тора с четырьмя степенями свободы. Исследование [10] базируется на допущении, что реакции подшипника — 
это квазилинейные силы с кубической нелинейностью. С учетом этого рассматривается воздействие радиально- 
го зазора на движение в пространстве динамически неуравновешенного ротора под влиянием сил внутреннего 
трения. В [11] показаны взаимосвязи поперечных и крутильных колебаний, возникающие при вращении ротора 
центрифуги. Разработана линеаризованная математическая модель ротора в упругих опорах, учитывающая вли- 
яние поперечных и крутильных колебаний. 

Проведенный обзор литературы позволяет утверждать, что исследование колебательного процесса и 
соответствующих параметров линейных механических систем представляет интерес. В данной работе рассмот- 
рен ротор, в частности его опорные узлы с зазорами. Цель исследования — определить значения под нагрузкой 
линейных эквивалентных жесткостных параметров названных узлов. 

Материалы и методы. В качестве ротора использована динамическая модель, предусматривающая 


вращение в упругих опорах с зазорами (рис. 1). 
Уа 
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Рис. |. Динамическая модель ротора, который вращается в упругих опорах с зазорами 
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На рис. 1 приняты следующие обозначения: | — опора ротора; 2 — граница контакта опоры с цапфой 
ротора; 3 — кривая относительного движения центра цапфы ротора; т — масса ротора; В (4) — радиальный 
зазор в роторных опорах (представлен функцией угла смещения его цапф от вертикальных направляющих); С,, 
С, — общие жесткости кожухов опорных узлов соответственно в горизонтальных и вертикальных направляю- 
ЩИХ; Уст 2. — статический сдвиг центра масс ротора в горизонтальном и вертикальном направлении, вызван- 
ный деформациями в опорных узлах; о — динамический угол наклона цапф ротора; о. — угол наклона от по- 
ложения равновесия, вызванный технологической нагрузкой на ротор; уч, 2. — полные динамические смещения 
центра масс ротора в указанных выше направлениях; / = Е2 / Ру — взаимосвязь вертикальной и горизонтальной 
компонент технологической нагрузки. 

Для определения потенциальной ( /Г) и кинетической (7) энергий исследуемой системы воспользуем- 
ся следующими равенствами: 


п-5С, |, +У, -6(а+а„)зт(а+а„)+5(о„ )зта | -- 


+5С. [2 +2„-6(“+а„)со$(о- 0 „)-+9(0„ )с0з а Гы 


—те | а +5 (а )соз и |; 
1 ее 1 
р = 5т(5. +2, +. 46 : 


Здесь А — полярный момент инерции, « — угловая скорость вращающегося ротора. 
Следует учитывать, что действующие силы не являются потенциальными. При этом обобщенные силы, 
которым соответствуют введенные ранее координаты у, 2, 9, примут вид: 


ЕЕ ==, 


а а 


[5(а-+а.„)соз(а+9„)+5'(а + ) т (а -+ а ) |+ 


О. =Е 


| +/[5(а+4„ эт (о „)-5' (а +, „ со (4+ 9 ) 


2 





Здесь 5' (© + о — производная радиального зазора по углу 9; } = — взаимосвязь вертикальной и гори- 


У 
зонтальной компонент технологической нагрузки. 
При данных условиях уравнения, выражающие динамику системы, можно записать: 


ту, +С, у, -С,6 (аа )зт(а+а)+С,6(а„ та „ = 0, 
т? +С.-2,-С.6(а-+а„ )соз(«-+а„)+С.6(о„ )соза „ = 0, 
т=|5(а+а„)зт(а +9 )-5' (+4 „)соз (+ 9 ) |- 
-С,| у, -5(а+а)зт(а+а„)+5 (0) то, |х 
х| 5'(а+а„)зт(а+а„)+5(а+0..„)со$ (0+9 ) |+ (1) 
+С. [2 -б(а-+о„)соз(а+а„)+5(а.„)созо.. |х 
х| 5 (+ )зт(а-+0„)- 6’ (а. )соз(а+ а) |= 
[5(а-+а.„)соз(а+9„)+5'(а+ „эт (+0 ) |+ 
° + 5(а+а„ т (о )-5' (о +0 „ )соз (4+0. | 
Рассмотрев более детально третье уравнение системы (1), определим угол наклона ©„, считая @&=0. 


После преобразований в уравнениях (1) получим более простую систему: 


ту, +2. у, + СХ, -2, 50, в 
ТР, + Си "2. + С. У. =0. 

Анализ уравнений (2) позволяет в первой аппроксимации значений рассматривать упругую опору с за- 
зором как опору с линейными упругими характеристиками (и в горизонтальных, и в вертикальных направляю- 
щих): 

а 


ЗИ У у __ 2 2у 
К, =С. о ВО "Ра Е. = Св "а Св "Уз: 


Здесь 


Машиностроение и машиноведение 


С^ 
> 


БИр://уезшК.Аоп$а.га 


Вестник Донского государственного технического университета. 2020. Т. 20, № 1. С. 61-67. 155М 1992-5950 е155М 1992-6006 
Тези! о{ оп ие Тесйтса! Оптуегзйу. 2020. Ио. 20, по. 1, рр. 61-67. 155М 1992-5980 е155М 1992-6006 





С. | те-г+ 03° м, (0+Х:2:2°.С. С.| те-г+с08° ©, (0+ Г.2.0°.С 
У У 


р ЕЕ НН 
ТА ТА 
С —С® = С.С. с05° о -0.2(0+1-2) 
экв экв 7 2 


и = 5" (а „)+25’(0„ ) —5(„)-5” (©. ) 2 
0=56(о„)-+5’(о„)- а, 2=0(9„)- 20 206) 
А=те-г+ сз а„(0+/.2)-(С..2*+С,.6?). (3) 


С", С” выражают эквивалентные характе- 


экв ? экв 


Результаты исследования. В системе (3) величины С” 


экв ? 
ристики жесткости опоры с зазором. 
Чтобы получить зависимости эквивалентных характеристик жесткости нелинейной механической си- 
стемы, воспользуемся разложением тригонометрических функций в системе (1), учитывая слагаемые более вы- 
сокого порядка. Для определения необходимых величин укажем: 


у. =А, Ш, 2, = В, -С05 Е, 
где В., А, — амплитуда общего колебания ротора в вертикальных и горизонтальных направляющих соответ- 


ственно, &, — частота колебаний. 


Таким образом, 
2п 


© “ | 
СУ т Ф, (у,,2, ) т ФЕ, 


экв 


У 
р (4) 
2 © р 
С, =—— | Ф. (у,,2, )с0з 8. 
пВ. 0 
Здесь Ф, (у,,2,), Ф. (у,,2,) — функции координат у,, 2, . Эти величины можно определить, обратившись к 


системе дифференциальных уравнений (1). 
Решая систему уравнений (4) методами интегрирования, представим величины в виде: 


С, (те"+С. с05 0, 5.2”) С’ 0 ди-0-5° (С? - 07.4, +С, -2° В?) 


ОЕ ВЕ ВЫВЕЗЕНЫ ЗЕ ВВНЫН -ННЕ ЗВННЫ 1 
А А (5) 

ве С, (тет + С, с0$ 0, 5.0) й С? сз „у. 2: (с; - 07. А, +С, -27 В?) 

” д 81 | 


Здесь 
Ь = с03° „(0+ 7-2), 
у =5(о.„)-+35'(а..„ во. —35"(о„)-5” (0%) 180, 
&=6(„ ) 120 „-36'(а„)-36" (0 „ ео „ +6" (0%). 
Анализируя (5), можно говорить о взаимозависимости не только между жесткостями корпусов в опорах 
ротора С,, С., но и амплитуд А,, В, его общих колебаний. 
Необходимо подчеркнуть, что жесткость С” выражается как взаимосвязь движений ротора в горизон- 


тальной и вертикальной плоскостях. В случае С>” =0 в системе (2) уравнения не будут связаны, что позволяет 
учесть все возможные варианты. Они перечислены ниже. 
Го —_ 7 — <“ 
— При отсутствии изменений радиального зазора со временем и © „=0 имеем С.’ =С.. Такой вари- 
ант характерен для холостого режима работы ротора. При этом радиальное смещение ротора по отношению к 


опоре в вертикальных направляющих является совсем небольшой величиной, хотя и более высокого порядка, 
чем смещение в горизонтальных направляющих. 


— Если принять 2=0,то С. =С.. Данная ситуация свойственна радиальному перемещению ротора 


по отношению к опоре в окружении точки, через которую проходит горизонтальная касательная к траектории 
относительного движения центра цапфы. 


д = 2 лее 
— При нулевом значении величины радиального зазора в опорах имеем С. =С, и С. =С.. 
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—ь о Е 
— Если отсутствуют изменения радиального зазора и ©„ = 90°, то С», =С,. Это возможно, если ро- 


тор работает при технологической нагрузке с подобранным зазором, т. е. в горизонтальных направляющих ра- 
диальное смещение принимает малые значения по сравнению с вертикальным направлением. 

— При 09=0 С.,, =С,. В данном случае радиальные смещения ротора имеют место в окрестности 
точки, через которую проходит вертикальная касательная к линии относительного движения центра цапфы. 

— Возможно выполнение равенства 9+ /.2=0. Это справедливо в случае радиального смещения ро- 
тора по отношению к опоре в окружении такой точки, через которую проведена нормаль, совпадающая с лини- 


ей действия результирующей сил Р, и Ё,. В этом случае С. =С, и С.„, =С.. 


Рассмотрим характеристики ротора на упругих опорах, жесткости которых равны С” 


экв ? 


2 
С*.. Обозна- 
чим = — эксцентриситет ротора, © — его угловая скорость. Уравнения колебаний представим в виде системы: 


.. у __ ря 2. 
ту +С- у, =т-5: 0 91101 


тЁ +С* -2 =Т:8- 67 СОЗФЁ (6) 
а экв ба — 
Решая (6), запишем выражения для амплитуд колебаний: 
т-=: 0” т-5: 0’ 
=, В= о. (7) 


у _ В оон 8 й ай 
С„-т-о С. -т- 0 


Э 


Заметим, что слагаемые С.„,, С. в знаменателях зависят от амплитуд А,, В, . Представим графически 


экв ? экв 


рассматриваемую систему (рис. 2). 


с’, (7, (90! 





В 


Рис. 2. Графическое представление динамических характеристик абсолютных колебаний ротора в опорах с зазорами 
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Поверхности Р., Р 


ЕЕ Р,. (резонансные) и С,, С. (скелетные) построены в предположении, что 


2 2 
амплитуды А, и В. изменяются независимо на основании (7) и равенств С., = т(®’ ) ‚С. = т(о: ) 
Амплитуды вынужденных колебаний А,, В, соответствуют частоте возмущающих сил ®,. Принимая 


во внимание этот факт, для получения выражений указанных амплитуд сначала построим плоскость ©’ — р : 


и Р 


122 


На линиях пересечения построенной плоскости с поверхностями Р., Р 


о Р,‚. будут точки с 


идентичными координатами. Эти точки — искомые амплитуды колебаний ротора. 


Л > “ 
Плоскости С, , С ” — скелетные поверхности линейной системы. В этом случае в выражениях для эк- 


2 


вивалентных жесткостей берутся только первые слагаемые. Они указывают на то, что частоты свободных коле- 
баний линейной системы не зависят от амплитуд. 
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Резонансные поверхности представляют колебания ротора в горизонтальных направляющих, а скелет- 
ные — в вертикальных. Они построены в различных координатных октантах. По форме скелетные поверхности 


являются эллиптическими параболоидами. На рис. 2 ©, ©. представляют собой частоты свободных колебаний 


системы. 

Обсуждение и заключения. Одной частоте могут подходить несколько точек. Это означает, что в ис- 
следуемой системе возможно несколько режимов колебаний, в том числе неустойчивых. Для перехода системы 
с одного устойчивого режима движения на другой необходимы некоторые внешние воздействия, что характер- 
но уже для нелинейных систем. 

В качестве примера применения полученных зависимостей можно привести задачу о вынужденных ко- 
лебаниях ротора вследствие его статической неуравновешенности. 

Представленные в работе зависимости (4) указывают на то, что эквивалентные угловые жесткости вза- 


имозависимы не только через жесткости С,, С, корпусов, но и через амплитуды А,, В, общих колебаний ро- 


тора. 

В результате исследования линейных характеристик опорных узлов нагруженного ротора получены 
следующие результаты. 

1. Предложена динамическая модель ротора, вращающегося в упругих опорах с зазорами. Данная мо- 
дель позволяет исследовать задачи определения линейных эквивалентных жесткостных характеристик опорных 
узлов. 

2. Проанализирована система уравнений и перечислены возможности применения формул. В частно- 
сти, с их помощью можно определить колебания ротора в горизонтальной плоскости на строгальных машинах, 
используемых для получения заготовок кожи. 
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